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Abstract 

We consider compensated spherical lens models and the caustic surfaces 
they create in the past light cone. Examination of cusp and crossover angles 
associated with particular source and lens redshifts gives explicit lensing models 
that confirm the claims of Paper I |||] , namely that area distances can differ by 
substantial factors from angular diameter distances even when averaged over 
large angular scales. 'Shrinking' in apparent sizes occurs, typically by a factor 
of 3 for a single spherical lens, on the scale of the cusp caused by the lens; 
summing over many lenses will still leave a residual effect. 



Subject headings: 
cosmology - gravitational lensing, -observational tests 



1 Introduction 



This paper continues the study started in Paper I |l| of how area distances behave in 
universes where strong gravitational lensing takes place. That paper considered the 
claim that although individual lensing masses alter area distances for ray bundles 
that pass near by them, photon conservation guarantees the same area distance- 
redshift relation as in exact^E It was shown in that this claimed compensation 
result is incorrect once one has passed caustics, which are necessarily the result of 
strong gravitational lensing; consequently (by continuity) the result is not true in 
general. Indeed it has to be wrong because area distances are determined by the 
gravitational field equations (essentially through the null Raychaudhuri equation) 
quite independently of the issue of photon conservation (which is determined by 
Maxwell's equations and is valid whatever the space-time curvature). Thus the latter 
cannot causally determine area distances. In fact at large distances, 'shrinking' takes 
place in that distant areas subtend smaller solid angles than they would in a FL 
universe model; and this effect will remain even when the observations relate to large 
angles. The way these small effects for individual lenses add up to give a significant 
averaged effect over the whole sky is discussed in Paper I. This may affect high-redshift 
number counts and Cosmic Background Radiation ('CBR') anisotropy observations 
at very small angles. 

The general argument has been given in I. Specific spherically symmetric examples 
(viewed from the centre, and so without caustics) are presented in [Q, which thereby 
gives a rigorous proof of the existence of the effect we claim; but the models used 
are unrealistic as models of the real universe. In this Paper (Paper II), we consider 
exact examples with caustics displaying the shrinking effects discussed in Paper I. We 
show how to calculate the magnitude of the effect analytically and numerically for 
single spherically symmetric compensated lenses to which we can apply the thin-lens 
approximation, and look in detail at 'top hat' lenses, which are the simplest in this 
class. We refer to and mainly follow the notation of Schneider, Ehlers and Falco 
(1992) [§ ('SEF'). 



2 The Gravitational Lensing Equations 

We consider compensated spherically symmetric lenses in an Einstein-de Sitter back- 
ground universe. The effect of the lens will be represented by the usual thin lens 
approximation, and we use the scaled variables of SEF. 



2.0.1 Background Relations 



The angular diameter distances between the observer and lens, observer and source, 
and lens and source in this background are D^, Dg and D^s respectively. In an 
Einstein - De Sitter model (f2 = 1, no clumping), 

„ _ 2c (1 + Zgyl^ (1 + Zs) - (1 + Zd) (1 + Zsf/^ 

Ho {l + z,){l + z,y 

[see SEF (4.57)], and Da (respectively, Dg) is obtainable from Dds by setting z^ — > 
0,Zs Zd (respectively, Zd 0). M follows, the dimensionless ratio TZ = ^ ^"^q"^' 
is important. For a given lens position Zd-, a.s Zg oc this has the limiting value 
"^oo = ^^\i+zj_)^'^ ' which has a maximum value of |^ = 0.21 when Zd = 9/7. 

If a source in a FLRW universe with scale factor a{t) emits a signal at time tg 
which is received at time t^, then the proper distance at time tg between the source 
and observer is i = a{tg) ^ . For an Einstein-de Sitter universe a{t) = t^/^, the 

Hubble parameter is H{t) = (2/3)t-\ and 1 + z = ^o^V^s^^ so 

i = 3ctl/^{tl^' - tY') = {2c/Ho){l + z)-^/^{VTT^ - 1) (2) 
is this distance in the background universe. 



2.1 Compensated lenses 

If the energy density is p{x) the fractional matter perturbation 5{x) in an inhomo- 
geneity is related to the matter source by 

S(x) = ^MH^ ^ p(f ) = po (5(f) + 1) . (3) 
Po 

where po is the average energy density over a hypersurface of constant time S, defined 
by 

p{x)d^x = Po . (4) 
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Is 

Integrating (|) over S, Jg 6{x)d^x = j- Is [p{x) - Po] d^x ; so by (|), 



po 

f3 



S{x)d^x = , (5) 

s 

which is the condition for a compensated perturbation that has been formed by 
rearrangement of matter in a uniformly distributed background with matter density 
Po- Equivalently, the density p averages out to the correct background value po] if 
this is not true, then the background density has been wrongly assigned |Q. Clearly 
this means that 6{x)IE positive. By (^), no negative densities will occur iff 

Po > , 6{x) > —1 everywhere. (6) 



We will assume these conditions to be true for the matter inhomogeneities causing 
lensing. Then in the lensing equations that follow, the quantities and relations will all 
refer to the variation from what they would have been in the background model (i.e. 
if there had been no lens). Thus the lensing 'surface mass density' cr will mean the 
projected surface mass density in the lens plane arising from the density difference 
6p = p — po from the background value, which will be chosen so that the compensation 
condition @ is true; and the 'bending angle' will be the deviation in direction at the 
lens from what it would have been in the background model. This will be given via 
the usual thin-lens equations, with the surface mass density as just defined. 



2.1.1 Simple compensated lenses 

We define a simple compensated lens ('SCL') to be a spherically symmetric compen- 
sated lens where 6{x) = S{\x\) is positive for an inner domain < < 1, negative 
for an outer domain 1 < |x| < A, (A > 1), and zero at larger radii, i.e. for A < 
This configuration would naturally arise by formation of a spherical massive object 
through gathering together material from an initially uniform substratum. In the 
sequel we consider a particularly simple form of SCL, namely a top hat lens. 



2.2 The lensing equation in scaled variables 

Given a choice of length scale .^o in the lens plane, there is a corresponding length 
scale r]Q = yt^o iii the source plane. From the position vector ff of the source relative 

to the optic axis in the source plane, and the impact vector ^ = D^O in the lens 
plane, where (vector) 9 is the observational angle from the optical axis, we define 
corresponding scaled variables x, y by 

^=f, y = ^ (7) 

[SEF (3.5)]. The surface mass density cr(|^|) for thin spherical lenses can be rescaled 
as 

a(^ox) c2 Hq 

«:(x) = ^, ^c.:=^^, x=\x\. (8) 

(SEF 5.4, 5.5). The lens equation can then be written in the very simple dimensionless 
form 

y = X — d{x) (9) 

(SEF 5.6, 8.6) where the scaled (vector) deflection angle a is related to the true 
(vector) deflection angle a by 

2cTZ - 

a{x) = ——a{^Qx) (10) 



(SEF 5.7). Because of the spherical symmetry, the deflection is radially inward and 
of magnitude y = \y\ given by 



X — aix) = X (llj 

X 



where 

m(x) 



2 r x'dx'nix') (12) 
Jo 



is the dimensionless mass m{x) within a circle of radius x (SEF 8.3); its first derivative 
is the dimensionless surface density (SEF 8.13): 

m' = = 2r/t(r) . (13) 

dr 

Also the change Ai in radial distance traveled by light in a given time, as measured 
at the source, is equal to the time delay caused by lensing [SEF (4.67), (5.45)][|, and 
can be rescaled to 

[SEF (5.11) and following]. This is given by 

Z = ^\a\^-ij{x), (15) 

where a is given by (|l^) and the (rescaled) deflection potential is 

ip(x) = 2 r x'K(x')ln ( -) dx' a(x) = . (16) 

Jo \x' J dx 

[SEF (8.7)-(8.9)]. 



2.2.1 The shrinking ratios 

Finally, pointwise over the sky, the angular shrinking factor 7 which relates observed 
distances corresponding to a given angle in the real lumpy universe to those in the 
background smoothed-out universe (see Paper I) is 

7 = \dyldx\ (17) 

which can be averaged over a stated angle AO to give the average angular shrinking 
factor (7) over that angle. Similarly the pointwise area shrinking factor f3 which 
relates observed areas corresponding to a given solid angle in the real lumpy universe 
to those in the background smoothed-out universe (see Paper I) is 

/3=|detJ|, J = 

This can be averaged over a solid angle Ai7 to give the average area shrinking factor 
{(3) over that sohd angle. 

"'^The time change calculated in these equations is at most a first order quantity, and so the 
change in radial distance travelled can be found from it to first order by calculating distance as if 
light travels on null geodesies in the background geometry. 



dy 
dx 



(18) 



2.2.2 The overall effect 



Together the radial and transverse equations ([T5|) , give the deflection of each null 
ray relative to the background geometry, and hence the shape of the perturbed light 
cone in the real lumpy space-time. These deflections are not independent: they are 
related by the fact that the speed of light is locally unity, so that the actual light 
path is stationary w.r.t. variation of the arrival time delay. Consequently a sideways 
deflection (which increases the distance to be traveled) is compensated by an inwards 
deflection (reducing the distance to be traveled), so the (tangential) lensing equation 
is a consequence of the (radial) time delay equation (SEF pp. 146-147, 170-171). It 
is this combination of radial and tangential effects, implied by the above equations, 
that gives the light cone caustics their characteristic shapes (see Figure 1). 



3 Angles and Distances 



Consider now angles and distances in the perturbed space-time. We start with an- 
gular diameter distance. Consider the source plane of a lens in direction (6). Image 
points of nearby directions will be displaced from their background position {i, x) by 
the (scaled) displacement (A£, Aff) {Z, y — x) = {Z, —a) where the flrst part is the 
radial component of the displacement, given by (p^D, and the second is the tangential 
component (in the source plane) given by ([TI|). If we move our viewing direction 
through an arc in the sky, the image point will move; for simplicity we will consider 
an arc where only one angular component 6 only varies. This gives a 2-dimensional 
section of the full S-dimensional^E 

As we vary 9 through d9, the (scaled) tangential distance traveled will be dy'^, and 
the radial change of distance dZ"^ will be much less than this. Thus the total distance 
traveled DT due to an angle increase AO is, to good approximation, 

dO , (19) 

where we sum all distances with a positive sign, thus determining the total increment 
in \y\ (see Paper 1). In terms of normalised magnitudes when spherical lensing takes 
place, the integrand is just \dy/d9\ = \1 — da / dx\dx / 06 . By contrast, the background 
distance is the same expression but with integrand \dx/d6\, and Distance gained DG 
is the distance moved from the starting point: 

DG{Ae) = [ = y{e) - y(0) . (20) 

In this case we subtract off those regions where dy/dO is negative, ending up simply 
with the increment of y. 



DT{A9) 



dm 

A9 89 



Now when is positive, distance traveled is the same as distance gained. How- 
ever when changes sign, we have cusps forming (see Paper 1) and in the formula 
( p!9D for distance traveled the integral is over the curve corresponding to all values of 
9 and hence traverses the cusp backwards and forwards, see Figure 1. This is different 
from distance gained; the latter is then given by the integral (|T9|) but where now the 
integral is over the (connected) curve 7 excluding the cusp sections, so that is 
positive over all the curve traversed. 

The Change in Distance Gained due to the presence of the deflector is small in 
all cases. The effect for large angular scales does not average to zero when we have a 
distribution of lenses, but it is very small (the change from the background value is 
given by the difference of —Dds « (^) at the two ends, corresponding to a minute of 
arc at most). The Change in Distance Traveled ADT is given by the integral (|19|), 
but now taken over all the closed loops 7c that are excluded when one calculates 
distance gained. The effect at each lens is small, but it is cumulative. Hence when 
there are a large number of lenses, the effect can be large (as discussed in Paper I). 

We are also interested in the true Cosmological Area Distance and soiE the ques- 
tion, as we look over a given solid angle, what area does that cover at the source? 
This change is given pointwise by the determinant of the lens equation, see ([T8|) . We 
give explicit expressions for this determinant in the following paragraphs. A radial 
increase of size will be partly compensated by a transverse decrease of size, see e.f. 
Gunn and Press 0), so the area distance will not relate very simply to the (radial) 
angular size distance. 

We see then that before caustics form, distance traveled and distance gained are 
both very similar (and very close to the background value, on large angular scales). 
Thereafter, they can be very different (as was argued in Paper I). To calculate this, 
we must locate the cusps and caustics. 

3.1 Caustics and Critical Curves 

Caustics in a source plane are the points in the plane where the Jacobean of the 
lensing map is singular. Critical curves are the points in the lensing plane where the 
light rays pass that will end up at caustics at the source plane. They can be located 
by determining the zeros of the Jacobean of the lensing map. The set of caustic points 
in space-time for all source planes form the space-time caustic set. 

3.1.1 The Jacobean 

Considering a spherical lens L centered at the origin of the Cartesian XF-plane, the 
Jacobian matrix J = ( ^ ] of the (transverse) lens mapping in the lens plane has 



determinant 

(SEF 8.16) which vanishes where either the first or the second brackets on the r.h.s. 
vanishes. 

When the first bracket in f^Tl) vanishes, the radius x is Xr such that 



= 1 ^ = 1 ^ = — . (22) 

Xj Xc Tc iZ 

Such a critical point occurs for example at r = {xc, 0); then there is a tangent vector 
St = (0, 1) to the critical curve at this point, and since the curve is tangential, Sj is an 
eigenvector with zero eigenvalue. Since the tangential critical curves are mapped onto 
the point r/ = in the source plane, there exists a caustic there which degenerates 
to a single point. The equation of the tangential critical curve in two dimensions in 
the lens plane is then simply x^ = x'^, where Xc solves (P^. This corresponds to an 
Einstein ring [many images, in a circle, of one point in the source plane]. 



The determinant det J in ([21| ) also vanishes where the lastiE 

da{xd) ^ 



(23) 



dx 

This equation describes radial critical curves. Again it corresponds to a circle in 
the lens plane. It has a radial eigenvector with eigenvalue zero. For instance, at 
r = (^,0), Er = (1,0). We see in the next section it corresponds to a caustic in the 
source plane [and a cusp in the surface of constant distance]. 



3.1.2 The Cross-over and Cusp Angles 

The lensing equation (^) is a two-dimensional vector equation with (transverse) com- 
ponents Hi and ?/2 while the radial equation gives the 3rd-component for the 2-d 
section of the past null cone in any surface of constant time. Lensing is radially in- 
ward with radial displacement magnitude given by ([TT|). The first term on the right 
is the position that would have been with no lens; the second term is the effect of the 
lens. 

To express this in terms of the observational angle 6 from the optical axis, we 
note from the relations ^ = D^O, x = ^/^q that x = DdO/^Q. Hence the magnitude 
equation takes the form 

y{e) = ^{^-^^ (^'^^)) • (24) 

(SEP 4.47b, 5.34). The cusp angles di and are determined by 



where again by the symmetry, 6i = —9^i. Differentiating (0), this occurs when 

= D,-Dds-^{Dde,)Dd (26) 

that is 

= i (27) 

determines the cusp angle di. These angles correspond to the radial critical points in 
equation (0). In terms of the bending angle diagram (Figure 2), this occurs where the 
curves y — x = yi are tangent to the curve a{x). The cusp physical size is DgOi, twice 
this distance is the difference between distance gained and distance traveled, to good 
approximation. JE 62 and 6^2 are related by 

y{e2) = y{e^2) = 0, (28) 

where by the spherical symmetry 62 = —9-2 and the self-intersection of the light cone 
(given by the first equality in this equation) occurs on the central line through the 
lens (as implied by the second equality). Thus we have from 



0-2 = ^&{D,e_2) (29) 

determines the cross-over angle 9_2- Thus the cross-over angles 62 and 6-2 correspond 
to the critical points satisfying equation (^2]). In terms of the bending angle diagram 
(see Figure 2) this occurs where the line y — x = Q intersects the curve a{x). 

An angle ^3 and corresponding impact parameter 0:3 yields the same image position 
as the cusp angle, on the other side of the caustic: y{x^) = —y{xi), and it is this 
angle that we treat as the cut-off in the caustic size. Henceforth, we refer to this as 
the cut-off angle 63 (and the cut-off on the other side occurs at 6-3 = —6*3). 

Finally, the maximum deflection caused by the lens occurs when 6 = ±6m, where 

^{D^e^) = (30) 

This does not correspond to either of the other angles; indeed it lies between them. 
For a SCL centred at ^ = 0, if cusps and cross-overs occur then generically 

< 01 < < 02 < ^3 . 



The two-dimensional picture obtained by suppressing one angular coordinate is 
as shown in Figure 1 (with one radial coordinate and one angular coordinate). Going 
to the full 3-dimensional picture, at the source plane, the whole picture is circularly 
symmetric about the optical axis at 6* = 0. The cross-over angles at = ±62 corre- 
spond to a circle in the lens plane but a point (a degenerate caustic) in the source 
plane; the cusp angles 6 = ±61 correspond to circles in both planes. 



We now apply the preceding theory to Top Hat models. 



4 The Top-Hat Matter Distribution 



density 5+ for < |x| < 1 and a constant outer density 5_ for 1 < |x| < A with A > 1. 
Then the compensation condition is 

6_ = -{\^-l)-H^, (31) 

Unless otherwise stated, we will assume that 5+ is positive (so 5_ is negative). Then 
the positivity condition (|) demands that 

< (5+ < (A^ - 1) ^ > (5_ > -1 , (32) 

using the scaled variables, and k{x) will take the form 

k{x)= C (A^v/T^^ - VA2 - x2 ) , 0<a;<l, (33) 

k{x)= -CVA2 - x2 , l<x<A, (34) 

k(x) = , X > A (35) 

where C = —2po6^/acr, with a central value k(0) = CA(A^ — 1) > and a junction 
value of /t(l) = — CVA^ — 1 < 0. The surface density will positive for < x < x+ < 1, 
negative for a;+ < a; < A, and zero for A < x, where 

x+ = \^{\'-l)/i\'-l) <1. (36) 



Substituting into (H) and integrating (p!2D to find the mass function m(x), we 
obtain the following: 

m{x)=Af{x), A = -^sJ-^, (37) 

O O'er 

where the function f{x) is given by 

/(x)= [l-^f-(l-x'f\ 0<x<l, (38) 

f{x)= l<x<X, (39) 

f{x)= 0, A<x. (40) 

The function /(x) is a continuous positive even function, with /(O) = df /dx{0) = 0, 
a single maximum value of f{xm) = (A^ — l)^/^/(A^ — 1)^^^ at x^ < 1 given by 
xl = A2(A^ - 1)/(A6 - 1) = xl, and junction values /(I) = (A^ - 1)3/2/A3, /(A) = 
= df/dx{X). Near zero it has the form 



It follows that m(x) is a continuous non-negative function with m(0) = and junction 
values m(A) = and m(l) = A{X'^ — l)^/^/A^. Its maximum value is at a; = x^, = x+, 
where m{xm) = Af{xm)- 

Consequently, because any SCL lens can be built up by a superposition of a 
sufficient number of top hat lenses, we see that the effective surface deflection mass 
M(r) is always positive and is exactly zero at the outer edge of the compensating 
region, that is the effective 2- dimensional surface density a is exactly compensated 
if the 3-dimensional fractional density 6 is precisely compensated. Hence there is no 
long-range effect due to the lens: precisely because it is correctly compensated, the 
deflection angle a = for impact parameters that lie outside x = X (where the density 
takes exactly the background value). Thus we note, (1) for compensated lenses, 
lensing effects occur only for rays that traverse the lens itself and its compensating 
region; (2) despite the negative values for a at some radii in such a compensated lens, 
the deflection angle is always positive. 

Collecting formulae resulting from (|1T],|12|) and (p7|-^), we have that for a spheri- 
cally symmetric top-hat matter distribution, 

a{x) = A^, (42) 

X 



where the constant is 



From (^) or (^), cross-overs occur where 

B{x) := ^ = 1 (44) 

X 

and from (|23| ) or (|27| ) caustics occur where 

da{x)/dx = 1, (45) 
The maximum bending angle am occurs where da/dx = 0. 

Consequently, 

(1) the bending angle a{x) is a continuous positive odd function with a{0) = 0, 
da/dx{0) = {3A/2) (^p^), a single maximum value at x^ < ■Co where xf^ = 
^S"-!)^ , and junction values a(l) = ^(A^ - l)3/2/A^ a(A) = = da/dx{X). 

(2) its slope da{x) / dx is an even continuous function with maximum value da{0) / dx 
(3y4/2A^)(A^ — 1) at the centre, positive from x = to Xm, negative from a; = to A, 
and zero thereafter, with junction values da{l)/dx = — (A/A^)(A^ -|- 2)y/X^ — 1 (here 
it takes its minimum value and its derivative d/dx{da/dx) is discontinuous, diverging 
from the left but finite on the right) and da{X)/dx = 0. Hence caustics occur iff 

{3A/2X^){X^-l)=A,ru>l (46) 



(with a degenerate case when equahty occurs). If they occur, say at x = X2 , then 
< a;2 < 1 and X2 satisfies 



da{x2) A 



dx X? 



+ 2x^)vTr^-(l + ^)Jl-^ =1, (47) 



with the corresponding angle 02 given by 02 = X2io/ Dd- 

(3) The function B{x) = a{x)/x = Af{x)/x^ (with f\x) given by (p^-|iO|)) is 
even, positive, monotonic decreasing, and continuous, with a maximum value 5(0) = 
(3A/2A2)(A2 - 1) at the centre, and junction values B{1) = iA/X){\'^ - l)^/yX, 
B{X) = 0. Hence cross-overs also occur iff ( ^) is satisfied. They can occur for any 
value of X > X2 up to A. If they occur, say at x = Xi, then 

Af{xi)/xl = l. (48) 

with the corresponding angle 6i given by 6i = xi^q/ Dd-, where f{x) is given either by 
(|38D (if Xi < 1) or by (|39|) (if 1 < xi < A) (one cannot tell a priori in which range 
it will lie; one has to try to solve one, and if there is no solution, solve the other). 
Then 62{zd, Zg) is the angle determining how large a part of the sky is covered by the 
Einstein circle corresponding to the cross-over surface z = Zg for lenses at Zd (giving 
multiple images of the central point at Zg). How this scales with Zd (for given z^) 
depends on how po, C,o, A and 71 scale with Zd. 

(4) Pointwise over the image, the area shrinking factor is given by /? = \det J\ given 
by (|2T|) . This can be evaluated from the formulae given above. Using the expansion 

flF) one can evaluate this determinant near the centre-line 6 = 0; the result is 

det J = (1 - Acrit) + O(x^) (49) 



which is 1 near the lens (when A is small) and goes to —Acrit (see (^61) ) when A is 
large (the minus sign because images are reversed). 

We can determine a value for the lensing parameter Mq either by directly estimat- 
ing the quantities in the definition (^3]), or by estimating the maximum bending angle 
am for lenses considered. For example, if A = 2, the r.h.s. of ( ^2]) has a maximum 
value of 0.70A (when x = 0.87). From (^3|) with the bending angle relation (|4^) and 
angle scaling relation (|I0|) we see that then Mq is determined by the relation 

Mo X 1.14x0.7=^^. (50) 

M 



4.1 Results 



We have written a series of Truebasic programmes that compute all the relevant 
quantities for Tophat lenses, as functions of (i) the determining parameter A, (ii) the 



source redshift Zs for fixed lens redshift Zd-, (iii) the lens redshift Zd for fixed source 
redshift Zg. We have experimented with parameter values that correspond to observed 
gravitational lensing systems; some of the results are given in the following tables and 
in Figures 3 to 5. This area shrinking ratio is about 3 after cusps have occurred, for 
scales of about 3 times the cusp scale, corresponding to the cusp image point where 
the deflection is the same size as at the cusp. 



4.2 Galaxy clusters 



We present a table of results for parameters corresponding to four well-known galaxy 
clusters that cause gravitational lensing (note that we are not making detailed models 
of these objects; rather we are using their observed properties to determine reasonable 
parameter values in our SCL model). From the cluster Abell 2218 (see refs. 0, H) 
we have selected as images the arcs at redshifts Zs = 2.6 and 3.3 respectively, as a 
case study, where the brackets imply this is evaluated at the angle cut-off angle. We 
then list the corresponding shrinking factor {(3) for these two images, at the cut-off 
angle, followed by their cusp and cross-over angles. The last column is the shrinking 
factor for the source placed at decoupling redshift Zg = 1200 . We also consider other 
lensing clusters Abell 963 (ref. [|| ), Abell 370 (refs. lll|, |2ll), and Abell 2390 (refs. 

mm)- 
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4.3 Galaxies 



We have also used a set of galactic lenses, as evidenced by multiple images of more 
distant objects, to provide parameters for our model, giving the second table. The 
first lens is often referred to as the 'clover leaf: +ifl413 + 117 has four images of a 
QSO at redshift Zd = 2.55 (See refs. [0, |2^). The second is the seen in QSO images 
A and B for the system 2345 + 007 correspond to a redshift Zd = 2.15, despite image 



A being 1.7mag brighter than image B (ref. |T^, The third is the triple radio 

source +MG201Q + 112 (see ref. The fourth is the QSO pair in 1635 + 267 with 

nearly equal redshift = 1.96 [|T3|. Finally a nearly full Einstein ring was observed 
in +MG1131 + 0456, albeit somewhat elliptic in shape (ref. ||^ ). 
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We find that the caustics shrinking factor tends to an average factor > 3 at large 
z (as required to get a 3-fold covering factor). However because of the divergence of 
the light rays within the cusps, it can be much larger for parameter values implying 
to strong lensing (the actual angle corresponding to the cusps is then much smaller 
than in the equivalent PL model). 



5 Conclusions 



Of particular interest is the way the cusp size and the "shrinking" vary with redshift of 
the source and of the lensing object. This depends on two things: firstly the variation 
of angular sizes with redshift, remembering (a) minimum angular apparent diameter 
occurs at z = 1.25, so that the maximum angle 9^ for cusps to form due to lenses of 
given size and strength will have minimum at that redshift; and (b) that the ratio of 
distances that enters cXcr saturates with increasing Zs (for given Zd) but has a maximum 
for each Zg at a Zd of about 0.6 which is thus the optimal distance for the lens in order 
to create cusps on the last scattering surface. M that are not typical of all galaxies 
or clusters; but they confirm in a concrete way the broad picture proposed in Paper 
I: an area 'shrinking' factor of 3 will occur for each lens that causes cusps, on the 
scale of the cusps (precisely: at the cut-off angle O^^ which gives the same deflection 



at the source as the cusp angle, but on the opposite side). The total effect when 
averaging over large angular scales will depend on what fraction of the sky is covered 
by these angles for all lenses at all smaller scales, as a function of redshift; some 
simple estimates of this overall effect were given in Paper I. To determine realistic 
multiplicity factors as a function of redshift will require simulations with multiple 
lensing and more realistic lens models, for example standard elliptical lens models 
determined by a velocity dispersion parameter and ellipticity parameters as in 
which allow an increase in the degree of multiple covering (because individual elliptic 
lenses can have a covering factor of 5). The effect will differ on angular scales, and 
will almost certainly be substantial due to micro-lensing, with an additional increase 
due to galactic and cluster lensing. The implication of this paper and Paper I is that 
it is incorrect to assume that areas average out to the background FL value on large 
angular scales; one can only know the true area ratios - expressed in the shrinking 
factors considered in this paper - by detailed calculation. 
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Figure 1: Shape of caustics in past light cone showing preferred geodesies and 
distance traveled. 




Figure 2: The bending angle diagram for two different redshifts. (a) one large, so 
other side of conjugate point Q; (b) One small, so this side of conjugate point Q. The 
number of images is the number of times the line y = x intersects the bending angle 
curve. Considering curve (a), firstly, there is one image corresponding to line (c); then 
there are two images for line (d) which is tangent to the curve and determines the 
cusp angle; there are three intersections for line (e) which determines the cross-over 
angles (as it corresponds to no displacement at the source plane); there are 3 images 
for generic position (f), again two images for line (g) as it passes through the cusp, 
and finally one image for line (h). Parameters based on the lens -l-MGl 131+0456. 




Figure 3: Variation in (a) cusp angle, (b) the crossover angle, and (c) the cut-off 
angle. The image redshift starts at the limiting value of Zs = = 0.231, and increases 
through the arc redshift of 0.914, up to the value Zg — 5. Parameters based on the 
lens A2390. 
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